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Analysis of drifter trajectories in the Gulf of Mexico has 
revealed the existence of a region on the southern portion of 
the West Florida Shelf (WFS) that is not visited by drifters 
that are released outside of the region. This so-called "for- 
bidden zone" (FZ) suggests the existence of a persistent 
cross-shelf transport barrier on the southern portion of the 
WFS. In this letter a year-long record of surface currents 
produced by a Hybrid-Coordinate Ocean Model simulation 
of the WFS is used to identify Lagrangian coherent struc- 
tures (LCSs), which reveal the presence of a robust and 
persistent cross-shelf transport barrier in approximately the 
same location as the boundary of the FZ. The location of 
the cross-shelf transport barrier undergoes a seasonal oscil- 
lation, being closer to the coast in the summer than in the 
winter. A month- long record of surface currents inferred 
from high-frequency (HF) radar measurements in a roughly 
60 km x 80 km region on the WFS off Tampa Bay is also 
used to identify LCSs, which reveal the presence of robust 
transient transport barriers. While the HF-radar-derived 
transport barriers cannot be unambiguously linked to the 
boundary of the FZ, this analysis does demonstrate the fea- 
sibility of monitoring transport barriers on the WFS using a 
HF-radar-based measurement system. The implications of 
a persistent cross-shelf transport barrier on the WFS for the 
development of harmful algal blooms on the shoreward side 
of the barrier are considered. 



1. Introduction 

Yang et al. [1999] presents the results of the analysis of 
trajectories of satellite-tracked drifters released during the 
period February 1996 through February 1997 on the conti- 
nental shelf in the northeastern portion of the Gulf of Mex- 
ico (GoM) . Inspection of the drifter trajectories depicted in 
Figure 2 of that paper reveals the presence of a trajectory- 
free triangular-shaped region on the southern portion of the 
West Florida Shelf (WFS), which has been referred to by 
those authors as a "forbidden zone" (FZ). Although little 
can be inferred about the spatio-temporal variability of the 
FZ from the aforementioned figure, the FZ almost certainly 
wobbles with a complicated spatio-temporal structure. This 
expectation finds some support in the seasonal analysis of 
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drifter trajectories presented by Morey et al. [2003], which 
also included trajectories of drifters released on the conti- 
nental shelf in the northwestern portion of the GoM. 

The presence of the FZ suggests the existence of a seem- 
ingly robust barrier on the WFS that inhibits the transport 
across the shelf. This cross-shelf transport barrier not only 
can have implications for pollutant dispersal, but may also 
be consequential for harmful algal blooms on the shoreward 
side of the barrier. 

In this letter we employ methods from dynamical systems 
theory to attain a twofold goal. First, we seek to demon- 
strate the robustness of the suggested cross-shelf transport 
barrier on the WFS in the inspection of drifter trajectories 
through the analysis of simulated surface currents. Second, 
we seek to demonstrate the feasibility of monitoring trans- 
port barriers on the WFS using high-frequency (HF) radar 
measurements. 

2. Lagrangian Coherent Structures (LCSs) 

Theoretical work on dynamical systems [e.g., Roller, 
2000; Holler and Yuan, 2000; Holier, 2001a, b, 2002; Shad- 
den et al., 2005] has characterized transport barriers in un- 
steady two-dimensional incompressible flows as Lagrangian 
coherent structures (LCSs). The theory behind LCSs will 
not be discussed in this letter. We note, however, that the 
LCSs of interest correspond to the stable and/or unstable 
invariant manifolds of hyperbolic points. 

An invariant manifold can be understood as a material 
curve of fluid, i.e., composed always of the same fluid par- 
ticles. Associated with a hyperbolic point in a steady flow 
there are two invariant manifolds, one stable and another 
one unstable. Along the stable (unstable) manifold, a fluid 
particle asymptotically approaches the hyperbolic point in 
forward (backward) time. Initially nearby fluid particle tra- 
jectories flanking a stable (unstable) manifold repel (attract) 
from each other at an exponential rate. These manifolds 
therefore constitute separatrices that divide the space into 
regions with dynamically distinct characteristics. Further- 
more, being material curves these separatrices cannot be 
traversed by fluid particles, i.e., they constitute true trans- 
port barriers. In an unsteady flow there are also hyperbolic 
points with associated stable and unstable manifolds. Un- 
like the steady case, these hyperbolic points are not still but 
rather undergo motion, which is typically aperiodic in pre- 
dominantly turbulent ocean flows. As in the steady case, the 
associated manifolds also constitute separatrices, and hence 
transport barriers, albeit in a spatially local sense and for 
sufficiently short time. The latter is seen in that generically 
there is chaotic motion in the vicinity of the points where 
the stable and unstable manifolds intersect one another after 
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Figure 1. Sequence of snapshots of FTLE field (1) computed backwards in time using surface currents generated by a 
HYCOM simulation of the WFS for year 2004. Maximizing ridges of FTLE field define attracting LCSs, which approximate 
unstable manifolds and thus constitute transport barriers. Note the presence of a triangular-shaped area on the southern 
portion of the WFS with small FTLEs bounded by the western Florida coast on the east, the lower Florida keys on the 
south, and large maximizing ridges of FTLE field on the west. The latter constitute a barrier for cross-shelf transport, 
which is seen to undergo a seasonal oscillation approximately about the western boundary of the FZ. 



successively stretching and folding. These intersections lead 
to the formation of regions with a highly intricate tangle-like 
structure. In these regions trajectories of initially nearby 
fluid particles rapidly diverge and fluid particles from other 
regions are injected in between, which constitutes a very 
effective mechanism for mixing. 

Identification of LCSs, which is impossible from naked- 
eye inspection of snapshots of simulated or measured veloc- 
ity fields and at most very difficult from the inspection of 
individual fluid particle trajectories, is critically important 
for understanding transport and mixing of tracers in the 
ocean. The computation of finite-time Lyapunov exponents 
(FTLEs) provides practical means for identifying repelling 
and attracting LCSs that approximate stable and unstable 
manifolds, respectively. The FTLE is the finite-time aver- 
age of the maximum expansion or contraction rate for pairs 
of passively advected fluid particles. More specifically, the 
FTLE is defined as 

^(xJ^irlnll&^xJII, (1) 

M 

where | denotes the L2 norm and </>t +T : x(t) 1— » x(i + r) 
is the flow map that takes fluid particles from their initial 
location at time t to their location at time t + r. The flow 



map <fil +r is obtained by solving the fluid particle motion 
regarded as a dynamical system obeying 

x = u(x,i), (2) 

where the overdot stands for time differentiation and u(x, i) 
is a prescribed velocity field. Repelling and attracting LCSs 
are then defined [Haller, 2002; Shadden et al, 2005] as max- 
imizing ridges of FTLE field computed forward (r > 0) and 
backward (r < 0) in time, respectively. 

We remark that while LCSs are fundamentally a La- 
grangian concept, the algorithm used to identify such struc- 
tures requires a high resolution Eulerian description of the 
flow. Recently, Lekien et al. [2005] has successfully applied 
this theory to identify transport barriers using HF-radar- 
derived surface velocity in the east Florida coast. 

3. LCSs Derived from Simulated Currents 

Numerical model output provides a flow description 
u(x,t) that is suitable for use to identify LCSs. Also, it 
has the advantage of allowing for a spatio-temporal coverage 
that is impossible to attain with existing observational sys- 
tems. Here we consider a year-long record of surface currents 
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produced by a Hybrid-Coordinate Ocean Model (HYCOM) 
simulation along the WFS for year 2004. 

The year-long record of simulated currents consists of 
daily surface velocity fields extracted in the WFS domain 
from a 0.04°-resolution, free-running HYCOM simulation 
of the GoM, itself nested within a 0.08°-resolution At- 
lantic basin data assimilative nowcast, which was gener- 
ated at the Naval Research Laboratory as part of a Na- 
tional Oceanographic Partnership Program in support of the 
Global Ocean Data Assimilation Experiment [Chassignet 
et al., 2006b, a]. The Atlantic nowcast was forced with re- 
alistic high-frequency forcing obtained from the U. S. Navy 
NOGAPS operational atmospheric model. It assimilated sea 
surface temperature and anomalous sea surface height from 
satellite altimetry with downward projection of anomalous 
temperature and salinity profiles. The nested GoM model 
was free-running and driven by the same high-frequency at- 
mospheric forcing. The topography used in both models 
was derived from the ETOP05 dataset, with the coastline 
in the GoM model following the 2 m isobath. Both models 
included river runoff. 

Figure 1 shows snapshots of FTLE field, which were com- 
puted using the software package MANGEN, a dynami- 
cal systems toolkit designed by F. Lekien that is available 
at http://www.lekien.com/~francois/software/mangen. At 
each time t the algorithm coded in MANGEN performs 
the following tasks. First, system (2) is integrated using 
a fourth-order Runge-Kutta-Fehlberg method for a grid of 
particles at time t to get their positions at time t + T, which 
are the values of the flow map at each point. This requires 
interpolating the velocity data, which is carried out employ- 
ing a cubic method. Second, the spatial gradient of the flow 
map is obtained at each point in the initial grid by central 
differencing with neighboring grid points. Third, the FTLE 
is computed at each point in the initial grid by evaluating 
(1). The previous three steps are repeated for a range of t 
values to produce a time series of FTLE field. Here we have 
set t — — 60 d so that the maximizing ridges of the FTLE 
fields shown in Figure 1 correspond to attracting LCSs. The 
choice r = —60 d was suggested by the time it takes a typical 
fluid particle to leave the WFS domain. Clearly, some par- 
ticles will exit the domain before 60 d of integration. In this 
case, MANGEN evaluates expression (1) using the position 
of each such particles prior exiting the domain. Note that 
due to the choice r = —60 d the time series of computed 
FTLE fields based on our year-long record of simulated cur- 
rents can only have a 10-month maximum duration. 

The regions of most intense red tonalities in each panel 
of Figure 1 roughly indicate maximizing ridges of FTLE 
field. These regions are seen to form smooth, albeit highly 
structured, curves that constitute the desired LCSs or trans- 
port barriers. Of particular interest is the triangular-shaped 
area on the southern portion of the WFS with small FTLEs 
bounded by the western Florida coast on the east, the lower 
Florida keys on the south, and large maximizing ridges of 
FTLE field on the west. The latter constitute a cross-shelf 
transport barrier that approximately coincides in position 
with the western boundary of the FZ identified in Yang et al. 



[1999]. This is most clearly visible during the period May 
through September. The sequence of snapshots of FTLE 
field in Figure 1 also reveals a seasonal movement of the 
cross-shelf transport barrier, being offshore during the win- 
ter and closer to the coast during the summer, which is in 
agreement with drifter observations [Morey et al., 2003]. 

4. LCSs Derived from Measured Currents 

Figure 2 shows a sequence of snapshots of FTLE field 
computed using surface velocity inferred from HF radar 
measurements taken during September 2003 in an approx- 
imately 60 km x 80 km domain on the WFS off Tampa 
Bay. 

The HF radar measurements consist of measurements 
made with Wellen radars, which mapped coastal ocean sur- 
face currents over the above domain with 30-minute inter- 
vals [Shay et al., 2006]. The computation of FTLEs was 
performed backward in time (r = —60 h) so that the maxi- 
mizing ridges of FTLE field in Figure 2 indicate attracting 
LCSs, which are analogous to perturbed unstable invariant 
manifolds. The numerical computation of the FTLEs was 
not carried out using the MANGEN software package. In- 
stead, we developed our own MATLAB codes, which, em- 
ploying a methodology similar to that outlined in the previ- 
ous section, allowed us to more easily handling FTLE com- 
putation based on velocity data defined on an irregular and 
totally open domain. 

The transport barrier character of the attracting LCSs 
identified in the FTLE fields depicted in Figure 2 is illus- 
trated by numerically simulating the evolution of three clus- 
ters of fluid particles. One of the clusters (black spot in the 
figure) is initially chosen on top of one LCS, while the other 
two clusters (dark-blue and magenta spots in the figure) are 
initially located at one side and the other of the LCS. Look- 
ing at the evolution of the clusters in time, we note that 
the black cluster remains on top of the LCS and stretches 
along it as time progresses. Also note that dark-blue and 
magenta clusters remain on two different sides, indicating 
the absence of flux across the LCS. 

We remark that LCSs are identifiable in the region cov- 
ered by the HF radar system through the whole month of 
September 2003. However, because of the limited domain of 
the radar footprint and the short deployment time window, 
we cannot say with certainty that any of the observed struc- 
tures corresponds to the boundary of the FZ. In spite of this 
uncertainty, our analysis of the HF radar measurements is 
highly encouraging inasmuch as it demonstrates the feasi- 
bility of tracking the evolution of the boundary of the FZ in 
near real time. 

5. Biological Implications 

In addition to being an interesting physical feature whose 
underlying dynamics deserves further study, the cross-shelf 
transport barrier on the WFS also has potentially important 
biological implications. 

The toxic dinofiagellate Karenia brevis has a rather re- 
stricted spatial distribution, primarily the GoM [Kusek 
et al., 1999]. While K. brevis exists in low concentrations in 
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Figure 2. As in Figure 1 but using HF-radar-derived surface currents in a roughly 60 km x 80 km domain on the WFS off 
Tampa Bay. The black, dark-blue, and magenta spots in each panel indicate clusters of passively advected fluid particles 
tracked from their initial location in the upper-left panel. Note how the black spot stretches along one of the identified 
LCSs and the dark-blue and magenta spots flanking this LCS are attracted to the LCS but do not cross it. 



vast areas of the GoM, it occasionally forms blooms along 
the northern and eastern coasts of the GoM [Wilson and 
Ray, 1956; Geesey and Tester, 1993; Tester and Steidinger, 
1997; Dortch et al, 1998; Kusek et at, 1999; Magana et al, 
2003]. The largest and most frequent blooms, however, oc- 
cur along the southern portion of the WFS [Steidinger and 
Haddad, 1981; Kusek et al., 1999]. Where these occur, there 
can be widespread death of fish, manatees, dolphins, turtles, 
and seabirds as a result of the brevetoxins produced by K. 
brevis [Bossart et al., 1998; Landsberg and Steidinger, 1998; 
Landsberg, 2002; Shumway et al., 2003; Flewelling et al., 
2005]. The brevetoxins also end up in an aerosol, which 
affects human respiration [Backer et at, 2003; Kirkpatrick 
et al., 2004]. Being slow growing while fast swimming algae 
[Chan, 1978; Brand and Guillard, 1981], dinoflagellates have 
the highest potential for achieving high population densities 
as a consequence of biophysical accumulation rather than 
actual growth. As a result of expatriate losses, K. brevis 
(and other dinoflagellates) would be expected to develop 
large blooms only in regions where stirring by currents is 
low. Indeed, many dinoflagellate blooms tend to occur in 
enclosed basins such as estuaries and lagoons, where expa- 
triate losses are reduced. However, as mentioned above, K. 
brevis often forms large blooms along the open coastline of 
the southern portion of the WFS, typically inside the FZ. 

We hypothesize that the cross-shelf transport barrier as- 
sociated with the FZ provides a mechanism that reduces K. 
brevis expatriation. A corollary of this hypothesis is that 
this barrier allows the nutrients from land runoff to accu- 



mulate near the coastline rather than being swept away by 
currents. While we cannot explain why K. brevis often dom- 
inates over other species of dinoflagellates in the FZ, we can 
predict that slow growing dinoflagellates will be more preva- 
lent within the FZ than outside. 

6. Summary 

In this letter we have shown that, when analyzed using 
dynamical systems methods, a year-long record of surface 
currents produced by a regional version of HYCOM reveals 
the presence of a cross-shelf transport barrier on the south- 
ern portion of the WFS, which is in approximately the same 
location as the boundary of the FZ identified earlier by Yang 
et al. [1999] using satellite-tracked drifter trajectories. The 
simulated cross-shelf transport barrier was robust, being 
present in all seasons while undergoing a seasonal oscilla- 
tion. The simulated cross-shelf transport barrier was closer 
to shore in the summer months than in the winter months 
in agreement with observations [Morey et al., 2003]. HF 
radar measurements were analyzed in a similar fashion and 
this analysis demonstrated the feasibility of experimentally 
monitoring transport barriers on the WFS using a system 
that can be operated in near real time. 
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